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ABSTRACT 
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other  fortifications  are  used  to  paralyze  the  forward  momentum  of  mechanized  anned 
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Army.  Due  to  the  Grizzly’s  mission  various  sensors;  laser  systems,  hydraulic  lines,  wires 
and  cameras  are  mounted  on  the  armor  hull,  which  are  exposed  to  various  types  of 
landmine  detonation  and  fragmentation.  This  thesis  studies  the  effects  of  shock  waves 
and  fragmentation  on  the  survivability  of  the  equipment  mounted  on  the  Grizzly’s 
armored  hull.  Models  of  an  OZM-72  antipersonnel  mine  are  developed  and  used  to 
simulate  the  detonation  and  fragmentation  phenomena.  The  analysis  results  obtained 
from  the  models  provide  a  basis  from  which  design  guidance  can  be  fonnulated  for 
protecting  equipment  or  personnel  from  this  threat. 
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I.  INTRODUCTION 


A.  BACKGROUND 

1.  Reason  for  Research 

Landmines  have  posed  a  longstanding  and  serious  threat  to  both  personnel  and 
vehicles  of  the  armed  forces.  By  using  complex  obstacles  such  as  mine  fields,  barbwire 
entanglements,  tank  ditches  and  other  fortifications  an  enemy  can  paralyze  the  forward 
movement  of  mechanized  armed  forces.  Due  to  these  threats  the  Grizzly  Breach  vehicle 
was  designed  for  the  United  States  Anny  [Ref.  8].  In  order  to  perfonn  its  mission  the 


Grizzly  Breach  Vehicle  has  various  sensors,  laser  systems,  hydraulic  lines,  wires  and 
cameras  mounted  on  its  armored  hull.  Since  they  are  out  side  the  protective  armor  hull 
the  various  sensors  and  equipment  are  exposed  and  vulnerable  to  the  detonation  and 
fragmentation  effects  of  bounding  landmines.  It  is  therefore  important  to  understand  the 
detonation  and  fragmentation  phenomena  of  bounding  landmines  in  order  to  provide 
guidance  in  the  protection  of  equipment  from  this  threat.  Modeling  and  simulation  of  this 
phenomena  is  conducted  using  an  OZM-72  bounding  antipersonnel  mine  to  provide 
insight  into  this  phenomena. 
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2.  Description  of  OZM-72  Antipersonnel  Mine  [Ref.  6] 

a.  General  Description: 

The  OZM-72  antipersonnel  mine  is  a  cylindrical  bounding  fragmentation 
mine  produced  by  the  former  Soviet  Union.  The  casing  of  the  mine  is  made  of  sheet 
steel,  which  contains  steel  wire  fragments  in  a  resin  matrix.  The  top  and  bottom  of  the 
mine  are  steel  plates  with  the  detonator  centrally  located  on  the  top  plate. 

b.  Method  of  Operation: 

The  mine  operates  by  first  igniting  a  propelling  charge,  which  lifts  the 
fragmentation  portion  of  the  mine  off  the  ground.  Once  the  mine  has  reached  a  height  of 
approximately  1.0  meter  an  anchor  wire  initiates  the  detonation  of  the  main  charge.  The 
main  charge  detonation  propels  the  steel  fragments  out  to  a  claimed  lethal  radius  of  25 
meters. 

c.  Technical  Specifications: 


Height  =160  mm 
Diameter  =  105  mm 
Total  Weight  =  5  kg 

Explosive  Weight  =  500  g  TNT 

Fragment  Dimensions: 

Length  =  5mm  Diameter  =  5mm 

Table  1.  OZM-72  Technical  Data 


Figure  2.  OZM-72  Landmine  [From  Ref.  5] 
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B. 


SCOPE  OF  REASEARCH 


This  paper  investigates  the  detonation  and  fragmentation  phenomena  of  bounding 
antipersonnel  mines.  This  is  undertaken  to  provide  a  better  understanding  of  the 
fragmentation  threat  posed  by  this  type  of  weapon.  To  accomplish  this  a  coupled 
Lagrangian  and  Eulerian  approach  was  used  in  modeling  this  phenomena.  Two- 
dimensional  and  three-dimensional  models  where  created  of  the  OZM-72  antipersonnel 
mine  using  the  AUTODYN  interactive  nonlinear  dynamics  analysis  software.  These 
models  were  used  in  providing  insight  into  the  detonation  and  fragmentation  phenomena. 
The  analysis  results  obtained  from  the  models  provide  a  basis  from  which  design 
guidance  can  be  formulated  for  protecting  equipment  or  personnel  from  this  threat. 
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II.  AUTODYN  COMPUTATIONAL  SETUP 


A.  SOLUTION  TECHNIQUE 

AUTODYN  solves  complex  dynamic  behavior  problem  by  discretization  of  the 
problem  in  both  time  and  space.  Time  is  broken  into  time  steps,  which  are  determined  by 
the  software  to  maintain  accuracy  and  stability  of  the  solution.  Space  is  divided  into 
elements,  which  are  six  sided  bricks  that  are  organized  by  a  structured  IJK  index  grid. 


Figure  3.  IJK  Brick  Element  [From  Ref.  4] 


Within  each  problem  several  IJK  index  grids  (sub  grids)  can  be  used  allowing  for 


complex  shapes  to  be  described  or  coupling  of  different  processor  types. 

B.  PROCESSOR  TYPE 

AUTODYN  has  several  processor  types  available,  which  can  be  used  in  a  given 
problem  to  provide  the  optimal  numerical  solution  for  the  given  problem.  The  processors 
used  in  this  study  are  the  Lagrange  and  the  Euler-FCT. 

1.  Lagrange 

The  Lagrange  processors  algorithms  are  based  on  the  finite  volume  method  used 
by  Wilkins  in  the  HEMP  code.  This  formulation  has  been  modified  slightly  in  order  to 
accommodate  forces  and  masses  at  the  nodes  similar  to  explicit  finite  element 
formulations.  This  modification  enabled  a  simpler  coupling  between  Lagrange  and  Euler 
processors  [Ref.  4].  The  Lagrange  processor  solves  problems  by  replacing  the  partial 
differential  equations  for  conservation  of  mass,  momentum  and  energy  with  finite 
difference  equations  solved  using  an  explicit  central  difference  scheme. 
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The  numerical  mesh  in  of  the  Lagrange  processor  moves  with  the  materials 
defined  in  the  sub  grid.  Meaning  that  no  material  is  transported  between  elements 
through  the  element  faces.  Properties  for  each  element  are  defined  either  at  the  node 
comers  or  at  the  cells  center,  as  shown  in  the  figure  below. 

>  x,u,F,m 
,p,e,p 


Figure  4.  Element  Property  Locations  [From  Ref.  3] 

Where  x  is  the  coordinate,  u  is  the  velocity,  F  is  the  force,  m  is  the  mass,  o  is  the  stress,  £ 
is  the  strain,  p  is  the  pressure,  e  is  the  internal  energy,  and  p  is  the  density. 

a.  Computational  Cycle 

The  steps  the  Lagrange  processor  conducts  for  each  time  step  in  a 
Lagrangian  sub  grid  are  shown  in  the  following  diagram. 


a.e 


Figure  5.  Lagrange  Computational  Cycle  [From  Ref.  3] 
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Since  the  algorithm  used  is  an  explicit  scheme  each  time  step  for 
Lagrangian  sub  grids  must  satisfy  the  CLF  or  Conart  condition 


At  < 


C 


(2.1) 


were  Ax  is  the  element  size  and  c  is  the  local  speed  of  sound.  The  minimum  value  for  the 
sub  grid  is  then  multiplied  by  a  safety  factor  of  two-thirds  and  is  used  to  advance  the 
solution  in  time.  The  reason  for  this  restriction  is  to  ensure  a  disturbance  dose  not 
propagate  across  an  element  in  a  single  time  step  and  to  maintain  stability  of  the 
algorithm. 

b.  Advantages 

The  advantages  of  using  the  Lagrange  processor  in  modeling  are  as 

follows: 


•  Computations  per  cycle  are  fewer  than  that  of  a  Eulerian  processor. 

•  Material  boundaries  and  interfaces  are  clearly  defined  and  do  not 
mix. 


•  Well  suited  for  modeling  solid  behavior  and  strength. 
c.  Disadvantages 

The  disadvantages  of  using  the  Lagrange  processor  in  modeling  are  as 

follows: 


•  Sever  element  distortions  lead  to  small  time  steps. 

•  Element  distortions  can  lead  to  grid  tangling  causing  the 
simulation  to  stop. 


These  disadvantages  can  be  overcome  by  the  use  of  the  erosion  feature  provided  in 
AUTODYN. 

2.  Euler  -  FCT 


The  Euler  -  FCT  processor  is  designed  specifically  to  solve  gas  dynamics 
problems  and  in  particular  blast  simulations.  FCT  stands  for  Flux  Corrected  Transport. 
The  algorithm  is  a  multi-dimensional  implementation  of  the  explicit  FCT  formulation  of 
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Boris  and  Book  [Ref.  4].  Like  in  the  Lagrange  processor,  the  Euler  processor  solves 
problems  by  replacing  the  partial  differential  equations  for  conservation  of  mass, 
momentum  and  energy  with  finite  difference  equations  solved  using  an  explicit  central 
difference  scheme. 

Unlike  the  Lagrange  processor,  the  Eulerian  numerical  grid  is  fixed  allowing 
materials  to  flow  through  the  element  faces  from  one  element  to  another.  Properties  for 
each  element  are  defined  either  at  the  node  corners  or  at  the  cells  center,  as  shown  in  the 
figure  below. 


.CUE 


,  x 

,p,e,p,u,m 


Figure  6.  Element  Property  Locations  [From  Ref.  3] 


Where  x  is  the  displacement,  u  is  the  velocity,  F  is  the  force,  m  is  the  mass,  o  is  the 
stress,  £  is  the  strain,  p  is  the  pressure,  e  is  the  internal  energy,  and  p  is  the  density. 

a.  Computational  Cycle 

The  steps  the  Euler  processor  conducts  for  each  time  step  in  a  Eulerian  sub 
grid  are  shown  in  the  following  diagram. 


Figure  7.  Eulerian  Computational  Cycle  [From  Ref.  3] 
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Since  the  algorithm  used  is  an  explicit  scheme  each  time  step  for  Eulerian 
sub  grids  must  satisfy  the  CLF  or  Conart  condition 


At  < 


(2.2) 


were  Ax  is  the  element  size,  c  is  the  local  speed  of  sound,  and  v  is  the  element  velocity. 
The  minimum  value  for  the  sub  grid  is  then  multiplied  by  a  safety  factor  of  two-thirds 
and  is  used  to  advance  the  solution  in  time.  The  reason  for  this  restriction  is  to  ensure  a 
disturbance  dose  not  propagate  across  an  element  in  a  single  time  step  and  to  maintain 
stability  of  the  algorithm. 

b.  Advantages 

The  advantages  of  using  the  Euler-FCT  processor  in  modeling  are  as 

follows: 


•  No  grid  distortions  or  tangling  reducing  time  step  or  stopping  the 
simulation. 

•  Accurate  higher  order  method  optimized  for  blast  type  problems. 
c.  Disadvantages 

The  disadvantages  of  using  the  Euler-FCT  processor  in  modeling  are  as 

follows: 


•  More  computations  per  cycle  are  required  than  that  of  a  Lagrange 
processor. 

•  Only  allowed  to  place  one  material  in  each  sub  grid. 

•  Only  allowed  to  use  the  ideal  gas  equation  of  state  to  define  the 
material  in  the  sub  grid. 


C.  MATERIAL  MODELS 

The  material  models  in  AUTODYN  can  be  broken  down  into  three  categories. 
These  categories  are  the  equation  of  state,  which  is  used  in  providing  volumetric  stress  or 
pressure.  The  strength  models,  which  are  the  constitutive  relations  for  detennining  the 
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deviatoric  stresses  by  Hooke’s  Law  and  a  plastic  yield  criteria.  The  last  is  the  failure 
model,  which  provides  a  criterion  for  determining  if  a  material  has  failed  and  no  longer 
has  strength. 

1.  Equations  of  State 

a.  Ideal  Gas 

The  ideal  gas  equation  of  state  relates  the  pressure  (p)  to  the  specific 
internal  energy  (e),  such  that 

p  =  (y-l)pe  (2.3) 

where  y  is  the  adiabatic  constant  (ratio  of  specific  heats)  and  p  is  the  density. 

b.  JWL 

The  JWL  equation  of  state  defines  the  explosive  shock  wave  pressure  (p ) 
as  a  function  of  the  specific  volume  (V)  and  the  specific  internal  energy  (e),  such  that 


p  -  C  |  1 - e 

rV 

V  rv  V 


+  C2  1--®-  + 

T  r2Vj  V 


where  variables  co,  Ci,  C2,  ri,  and  ro  are  constants  for  a  specific  explosive.  These 
constants  are  determined  from  data  collected  during  cylinder  test. 

c.  Linear 


The  linear  equation  of  state  considers  pressure  (p )  a  function  of  density  (p) 
by  an  approximation  to  Hook’s  Law,  such  that 


P  =  Kju 


where  p  =  \yp  J  - 1  the  compression,  p0  is  the  reference  density,  and  K  is  the  materials 
bulk  modulus. 


2.  Strength  Models 


The  main  purpose  of  the  strength  models  is  to  determine  when  a  material  should 
yield.  In  AUTODYN  all  the  strength  models  use  the  Von  Mises  yield  criterion  to 
determine  the  elastic  limit  and  the  onset  of  plastic  flow.  The  Von  Mises  yield  criteria 
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states  that  the  deviatoric  stresses  (distortion)  si,  S2,  and  S3  are  related  to  the  yield  strength 
(JO  by 

222  2Y2 

s\  +s2  +  S3  =  ^  (2.6) 

meaning  the  onset  of  yielding  is  purely  a  function  of  the  distortion  of  the  material  and 
does  not  depend  on  the  local  hydrostatic  pressure  unless  the  yield  strength  is  a  function  of 
pressure. 

a.  Von  Mises  Model 

The  Von  Mises  model  is  the  simplest  and  most  convenient  criteria  to 
describe  the  elastic  limit  and  the  transition  to  plastic  flow.  This  model  defines  the  yield 
stress  (Y)  as  a  constant,  which  is  entered  by  the  user. 

3.  Failure  Model 

a.  Bulk  Strain 

This  simple  model  states  that  bulk  failure  of  a  material  element  will  occur 
when  the  effective  plastic  strain  exceeds  a  specified  value.  Once  this  limit  is  reached  the 
stress  deviators  are  set  to  zero  and  the  element  material  can  no  longer  sustain  any  shear 
strength. 

4.  Erosion  Model 

During  simulations  some  Lagrangian  element  can  become  severely  distorted 
which  can  slow  the  progress  of  the  calculations.  In  AUTODYN  these  cells  can  be 
removed  from  the  calculation  if  a  predetermined  strain  exceeds  a  user  specified  limit. 

D.  PROCESSOR  COUPLING 

Due  to  AUTODYN’s  coupling  ability,  separate  numerical  sub  grids  with  the  same 
or  different  types  of  processors  can  be  used  to  describe  a  given  problem.  This  allows  for 
the  use  of  the  best  type  of  processor  to  be  used  in  modeling  the  problem.  Since  the 
interactions  of  the  sub  grids  are  also  need  to  be  calculated  the  following  is  a  diagram  of 
the  complete  computational  cycle. 


11 


Figure  8.  Complete  Computational  Cycle  [From  Ref.  7] 


For  this  study  the  types  of  processor  coupling  that  were  used  are  the  joined 
Lagrange,  impact/slide  interface  and  the  Euler-Lagrange  coupling. 

1.  Joined  Lagrange 

For  the  joined  Lagrange  type  of  processor  coupling  the  connected  nodes  of 
Lagrangian  sub  grids  are  fused  together  and  are  regarded  as  a  single  node  in  the 
calculations.  Meaning  the  joined  nodes  will  remain  together  unless  the  joined  condition 
is  removed  or  element  erosion  occurs. 

2.  Impact/Slide  Interface 

Impact/Slide  interface  is  used  in  determining  the  interactions  between  Lagrangian 
sub  grids  which  impact  or  slide  along  each  other.  This  interaction  algorithm  is  based  on 
the  use  of  a  small  gap  to  detennine  if  sub  grids  are  interacting.  This  gap  defines  a  zone 
around  each  interacting  face  that  is  used  in  detennining  if  the  sub  grids  will  collide 
during  a  time  step.  To  determine  if  sub  grids  interact  an  array  of  the  current  surface  faces 
of  one  sub  grid  is  comparing  to  the  surface  nodes  of  another  sub  grid.  If  any  of  the 
surface  nodes  enter  the  other  grids  gap  zone,  momentum-conserving  interactions  are 
computed  to  prevent  the  penetration  of  the  sub  grids.  This  procedure  is  repeated  by 
exchanging  the  sub  grids  in  order  to  provide  symmetry  to  the  process  [Ref.  3] 
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3.  Euler  -  Lagrange  Coupling 

The  Eulerian  and  Lagrangian  sub  grids  are  coupled  by  using  an  algorithm,  which 
defines  a  stress  profile  for  the  Lagrangian  sub  grid.  The  stress  profile  is  based  on  the 
geometric  flow  constraint  placed  on  the  Eulerian  sub  grid  as  the  Lagrangian  grid  passed 
through  it.  While  Lagrangian  sub  grid  moves  and  distorts  through  the  Eulerian  sub  grid 
elements  are  covered  and  uncovered.  Due  to  this  fact  a  sophisticated  logic  is  used  to 
avoid  the  Eulerian  elements  from  becoming  to  small.  This  is  done  in  order  to  ensure  the 
time  step  is  not  severely  reduced.  With  this  type  of  coupling  algorithm  complex  fluid- 
structure  interaction  problems  can  be  solved  using  the  best  processor  for  the  job. 

E.  REMAPPING 

Remapping  is  a  method  by  which  a  user  is  able  to  take  the  solution  of  one  analysis 
and  impose  it  upon  all  regions  of  a  selected  region  of  a  different  model.  This  is 
accomplished  by  creating  a  remap  data  fill,  which  contains  the  data  of  the  analysis  that  is 
to  be  remapped.  The  file  is  then  loaded  into  the  receptor  sub  grid,  which  provides  the 
receptor  sub  grid  with  the  data  from  the  previous  solution. 
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III.  VULNERABILITY  ASSESSMENT 


A.  DAMAGE  MECHANISM 


For  bounding  antipersonnel  mines  the  damage  is  created  by  the  detonation  and 
fragmentation  phenomena.  In  these  phenomena  the  energy  of  an  explosive  shock  wave  is 
transferred  to  the  fragmentation  material,  which  is  used  to  damage  a  target.  For 
antipersonnel  mines  the  explosive  shock  wave  is  produced  by  the  rapid  energy  release  of 
an  explosive  charge.  The  shock  wave  produced  by  this  energy  release  travels  faster  than 
the  speed  of  sound,  has  a  limited  duration  and  possesses  a  significant  amount  of  energy. 
The  energy  contained  within  the  shock  wave  is  then  transferred  to  the  fragmenting 
material  in  the  form  of  kinetic  energy.  The  kinetic  energy  of  the  fragmenting  material  is 
the  mechanism  by  which  damage  is  caused  to  targets.  The  kinetic  energy  of  a  fragment  is 
found  by  the  equation: 


(3.1) 


where  m  is  the  mass  of  the  fragment  and  v  is  the  velocity  of  the  fragment.  Since  the 
fragments  travel  thru  the  air  the  velocity  of  the  fragments  will  decrease  due  to  the  drag 
created  by  the  air.  The  velocity  of  a  fragment  at  a  specified  distance  can  be  determined 
by  the  equation  [Ref.  1 1] 


-pCdAS 

v0e  2m 


(3.2) 


where  v0  is  the  initial  velocity,  p  is  the  density  of  the  fluid,  Cd  is  the  drag  coefficient  of 
the  fragment,  A  is  the  cross-sectional  area  of  the  fragment,  S  is  the  distance,  and  m  is  the 
mass  of  the  fragment. 

B.  PROBABILITY  OF  KILL 


The  probability  of  kill  (/\)  is  a  statistical  measure  of  a  targets  ability  to  withstand 
the  damage  caused  by  one  or  more  damage  mechanisms.  Probability  of  kill  is  a  function 
of  whether  a  damage  mechanism  reaches  a  target  and  the  likelihood  the  target  will  be 
destroyed  or  killed  by  the  damage  mechanism.  In  this  study  the  probability  of  kill  was 
determined  based  on  the  damage  mechanism  described  in  section  A  of  this  chapter. 
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For  antipersonnel  mines  the  probability  of  kill  is  determined  by  finding  the 
number  of  fragments  that  impact  a  target  and  the  probability  of  kill  per  fragment  (Pyf). 
From  these  values  a  probability  of  kill  is  detennined  by  the  use  of  the  following 
equations  [Ref.  2  and  1 1] 

A  =  1  -  (l  —  Pw  )N‘  For:  N  >  1  (3.3) 

Pk  =  NXPm  )  For:  0  <  N  <  1  (3.4) 

where  N,  is  the  number  of  fragments  impacting  the  target. 

For  this  study  the  probability  of  kill  is  determined  by  finding  the  average  kinetic 
energy  of  the  fragments  from  the  average  velocity  of  the  fragments.  From  this  the 
probability  of  kill  per  fragment  was  determined  by  linear  interpolation  of  Table  2. 


Kinetic  Energy  (Joules) 

Target 

>-b 

II 

© 

Pk/f  =  0.5 

£ 

ll 

o 

to 

Soft 

100 

1000 

4000 

Hard 

4000 

10000 

20000 

Table  2.  Probability  of  Kill  Per  Fragment  [After  Ref.  1 1] 


Where  soft  targets  can  be  considered  to  be  personnel  or  fragile  sensors  and  hard  targets 
can  be  considered  to  be  body  armor  or  light  metals  like  aluminum. 


The  number  of  fragments  impacting  a  target  is  detennined  by  considering  an  ideal 
case.  In  this  case  the  mine  is  considered  to  be  level,  the  spray  pattern  of  the  fragments  is 
cylindrical  with  the  fragments  only  traveling  in  the  radial  direction.  The  targets  are 
considered  to  be  in  the  direct  path  of  the  fragments  and  have  a  height  equal  to  that  of  the 
mine.  The  following  equation  is  used  to  determine  the  number  of  fragments  impacting 
the  target 


N  =  w 


N 


2nR 


(3.5) 


where  w  is  the  width  of  the  target,  Nm  is  the  number  of  fragments  in  the  mine,  and  R  is  the 
distance  the  target  is  from  the  mine. 
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IV.  MINE  EXPLOSION  MODELS 


A.  EXPLOSIVE  MODEL 

The  explosive  for  the  mine  models  was  modeled  separately  from  the  remainder  of 
the  mine.  The  results  produced  by  the  explosive  model  were  placed  into  the  two- 
dimensional  and  three-dimensional  mine  models  using  the  remap  feature  available  in 
AUTODYN.  This  was  accomplished  in  order  to  reduce  the  computation  time  for  the 
model  and  so  that  the  shock  wave  pressure  could  be  modeled  using  the  JWL  equation  of 
state. 

The  explosive  model  is  composed  of  a  single  sub  grid,  which  used  the  Eulerian 
processor  and  was  filled  with  a  material  model  that  has  a  JWL  equation  of  state.  The 
material  model  input  deck  can  be  found  in  Appendix  A.  The  sub  grid  is  composed  of  a 
square  that  is  67mm  in  the  x  direction  and  27mm  in  the  y  direction  with  an  I  index  of  68 
and  a  J  index  of  28  (1876  elements).  In  addition,  an  axial  symmetry  boundary  and  a 
detonation  boundary  were  placed  along  the  x-axis.  Figure  9  is  a  depiction  of  the  TNT 
sub  grid. 


Y 


A 


Figure  9.  Explosive  Model 


The  dimensions  of  the  TNT  sub  grid  were  chosen  in  order  to  provide  the  appropriate 
mass  of  explosive  to  the  models. 
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B. 


TWO-DIMENTIONAL  MINE  MODELS 


The  two-dimensional  mine  models  were  developed  based  on  a  quarter  of  a  slice  of 
the  mine  when  viewed  from  the  top  or  bottom  of  the  mine.  Each  of  these  models 
contains  three  sub  grids  and  has  symmetry  boundaries  along  the  x,  y  and  z-axis. 

1.  Four-Row  Mine  Model 

Figure  10  is  a  depiction  of  the  assembled  four-row  mine  model. 


Figure  10.  Four-Row  Mine  Model 


In  Figure  10  the  blue  portion  is  the  air  sub  grid,  the  gray  is  the  shell  sub  grid  and  the 
green  is  the  wire  sub  grid. 

a.  Air  Sub  Grid 

The  air  sub  grid  used  the  Euler-FCT  processor  and  was  filled  with  a 
material  model,  which  used  the  ideal  gas  equation  of  state.  The  material  input  deck  can 
be  found  in  Appendix  B.  The  sub  grid  is  composed  of  a  block  that  is  5mm  in  the  x 
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direction  and  300mm  in  both  the  y  and  z  directions  with  an  I  index  of  6,  and  a  J  and  K 
index  of  101  (50000  elements).  In  addition  non-reflecting  boundaries  were  placed  along 
the  faces  indicated  in  Figure  1 1 . 

Non  -  Reflecting 


b  Z 

Figure  1 1.  Four-Row  Model  of  the  Air  Sub  Grid 


The  explosive  shock  wave  produce  by  the  explosive  model  was  remapped  into  the  air  sub 
grid  after  the  sub  grid  was  produced. 

b.  Wire  Sub  Grid 

The  wire  sub  grid  used  the  Lagrange  processor  and  was  fdled  with  a 
material  model,  which  used  the  linear  equation  of  state.  The  material  input  deck  can  be 
found  in  Appendix  C.  The  sub  grid  is  composed  of  34  cylinders  each  having  a  diameter 
of  5mm,  a  length  of  5mm  and  arranged  as  shown  in  Figure  12.  Each  of  the  cylinders  has 
an  I  and  J  index  of  5  and  a  K  index  of  4  (48  elements). 
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Figure  12.  Four-Row  Model  of  the  Wire  Sub  Grid 

c.  Shell  Sub  Grid 

The  shell  sub  grid  used  the  Lagrange  processor  and  was  fdled  with  a 
material  model,  which  used  the  linear  equation  of  state,  the  Von  Mises  strength  model, 
the  bulk  strain  failure  model  and  the  erosion  model.  The  material  input  deck  can  be 
found  in  Appendix  D.  The  sub  grid  is  composed  of  14  of  a  hollow  cylinder  with  an  inner 
radius  of  50mm,  an  outer  radius  of  52mm  and  a  length  of  5mm  with  a  I  index  of  2,  a  J 
index  of  51,  and  a  K  index  of  4  (153  elements).  Figure  13  is  a  depiction  of  the  shell  sub 
grid. 


Figure  13.  Four-Row  Model  of  the  Shell  Sub  Grid 
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2. 


Three-Row  Mine  Model 


Figure  14  is  a  depiction  of  the  assembled  three-row  mine  model. 


Figure  14.  Three-Row  Mine  Model 

In  Figure  14  the  blue  portion  is  the  air  sub  grid,  the  gray  is  the  shell  sub  grid  and  the 
green  is  the  wire  sub  grid. 

a.  Air  Sub  Grid 

The  air  sub  grid  used  the  Euler-FCT  processor  and  was  filled  with  a 
material  model,  which  used  the  ideal  gas  equation  of  state.  The  material  input  deck  can 
be  found  in  Appendix  B.  The  sub  grid  is  composed  of  a  block  that  is  5mm  in  the  x 
direction  and  90mm  in  both  the  y  and  z  directions  with  an  I  index  of  5  and  a  J  and  K 
index  of  46  (8100  elements).  In  addition  non-reflecting  boundaries  were  placed  along  the 
faces  indicated  in  Figure  15. 
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Figure  15.  Three-Row  Model  of  the  Air  Sub  Grid 


The  explosive  shock  wave  produce  by  the  explosive  model  was  remapped  into  the  air  sub 
grid  after  the  sub  grid  was  produced. 

b.  Wire  Sub  Grid 

The  wire  sub  grid  used  the  Lagrange  processor  and  was  fdled  with  a 
material  model,  which  used  the  linear  equation  of  state.  The  material  input  deck  can  be 
found  in  Appendix  C.  The  sub  grid  is  composed  of  30  cylinders  each  having  a  diameter 
of  5mm,  a  length  of  5mm  and  arranged  as  shown  in  Figure  16.  Each  of  the  cylinders  has 
an  I  and  J  index  of  5  and  a  K  index  of  2  (16  elements). 
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Figure  16.  Three-Row  Model  of  the  Wire  Sub  Grid 
c.  Shell  Sub  Grid 

The  shell  sub  grid  used  the  Lagrange  processor  and  was  fdled  with  a 
material  model,  which  used  the  linear  equation  of  state,  the  Von  Mises  strength  model, 
the  bulk  strain  failure  model  and  the  erosion  model.  The  material  input  deck  can  be 
found  in  Appendix  D.  The  sub  grid  is  composed  of  14  of  a  hollow  cylinder  with  an  inner 
radius  of  50mm,  an  outer  radius  of  52mm  and  a  length  of  5mm  with  a  I  index  of  2,  a  J 
index  of  21,  and  a  K  index  of  4  (60  elements).  Figure  17  is  a  depiction  of  the  shell  sub 
grid. 


Figure  17.  Three-Row  Model  of  the  Shell  Sub  Grid 
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3. 


Two-Row  Mine  Model 


Figure  18  is  a  depiction  of  the  assembled  two-row  mine  model. 


In  Figure  18  the  blue  portion  is  the  air  sub  grid,  the  gray  is  the  shell  sub  grid  and  the 
green  is  the  wire  sub  grid. 

a.  Air  Sub  Grid 

The  air  sub  grid  used  the  Euler-FCT  processor  and  was  filled  with  a 
material  model,  which  used  the  ideal  gas  equation  of  state.  The  material  input  deck  can 
be  found  in  Appendix  B.  The  sub  grid  is  composed  of  a  block  that  is  5mm  in  the  x 
direction  and  100mm  in  both  the  y  and  z  directions  with  an  I  index  of  3  and  a  J  and  K 
index  of  36  (2450  elements).  In  addition  non-reflecting  boundaries  were  placed  along  the 
faces  indicated  in  Figure  19. 
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Figure  19.  Two-Row  Model  of  the  Air  Sub  Grid 


The  explosive  shock  wave  produce  by  the  explosive  model  was  remapped  into  the  air  sub 
grid  after  the  sub  grid  was  produced. 

b.  Wire  Sub  Grid 

The  wire  sub  grid  used  the  Lagrange  processor  and  was  fdled  with  a 
material  model,  which  used  the  linear  equation  of  state.  The  material  input  deck  can  be 
found  in  Appendix  C.  The  sub  grid  is  composed  of  17  cylinders  each  having  a  diameter 
of  5mm,  a  length  of  5mm  and  arranged  as  shown  in  Figure  20.  Each  of  the  cylinders  has 
an  I  and  J  index  of  5  and  a  K  index  of  2  (16  elements). 
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Figure  20.  Two-Row  Model  of  the  Wire  Sub  Grid 

c.  Shell  Sub  Grid 

The  shell  sub  grid  used  the  Lagrange  processor  and  was  fdled  with  a 
material  model,  which  used  the  linear  equation  of  state,  the  Von  Mises  strength  model, 
the  bulk  strain  failure  model  and  the  erosion  model.  The  material  input  deck  can  be 
found  in  Appendix  D.  The  sub  grid  is  composed  of  14  of  a  hollow  cylinder  with  an  inner 
radius  of  5 1 .5mm,  an  outer  radius  of  52.5mm  and  a  length  of  5mm  with  a  I  index  of  2,  a  J 
index  of  1 1  and  a  K  index  of  2  (10  elements).  Figure  21  is  a  depiction  of  the  shell  sub 
grid. 


Figure  21.  Four-Row  Model  of  the  Shell  Sub  Grid 
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c. 


THREE-DIMENTIONAL  MINE  MODEL 


The  three-dimensional  mine  model  was  developed  based  on  one  eighth  of  the 
mine  and  had  a  wire  sub  grid  similar  to  that  of  the  two  row  mine  model.  Unlike  the  two- 
dimensional  mine  models  the  three-dimensional  mine  model  has  an  additional  sub  grid 
which  models  the  top  or  bottom  of  the  mine.  Like  the  two-dimensional  models  symmetry 
boundaries  were  created  along  the  x,  y,  and  z-axis.  Figure  22  is  a  depiction  of  the 
assembled  three-dimensional  mine  model. 


Figure  22.  Three-Dimensional  Mine  Model 

1.  Air  Sub  Grid 

The  air  sub  grid  used  the  Euler-FCT  processor  and  was  filled  with  a  material 
model,  which  used  the  ideal  gas  equation  of  state.  The  material  input  deck  can  be  found 
in  Appendix  B.  The  sub  grid  is  composed  of  a  block  that  is  1 10mm  in  the  x  direction  and 
90mm  in  both  the  y  and  z  directions  with  an  I  index  of  46  and  a  J  and  K  index  of  36 
(55125  elements).  In  addition  non-reflecting  boundaries  were  placed  along  the  block 
faces  indicated  in  Figure  23. 
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Non  -  Reflecting 
Boundary 


Figure  23.  Three-Dimensional  Model  of  the  Air  Sub  Grid 

The  explosive  shock  wave  produce  by  the  explosive  model  was  remapped  into  the  air  sub 
grid  after  the  sub  grid  was  produced. 

2.  Wire  Sub  Grid 

The  wire  sub  grid  used  the  Lagrange  processor  and  was  fdled  with  a 
material  model,  which  used  the  linear  equation  of  state.  The  material  input  deck  can  be 
found  in  Appendix  C.  The  sub  grid  is  composed  of  153  cylinders  each  having  a  diameter 
of  5mm,  a  length  of  5mm  and  arranged  as  shown  in  Figure  24.  Each  of  the  cylinders  has 
an  I  and  J  index  of  5  and  a  K  index  of  2  (16  elements). 
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Figure  24.  Three-Dimensional  Model  of  the  Wire  Sub  Grid 

3.  Shell  Sub  Grid 

The  shell  sub  grid  used  the  Lagrange  processor  and  was  filled  with  a  material 
model,  which  used  the  linear  equation  of  state,  the  Von  Mises  strength  model,  the  bulk 
strain  failure  model  and  the  erosion  model.  The  material  input  deck  can  be  found  in 
Appendix  D.  The  sub  grid  is  composed  of  %  of  a  hollow  cylinder  with  an  inner  radius  of 
5 1 .5mm,  an  outer  radius  of  52.5mm  and  a  length  of  70mm  with  a  I  index  of  2  and  a  J  and 
K  index  of  1 1  (100  elements).  Figure  25  is  a  depiction  of  the  shell  sub  grid. 


Figure  25.  Three-Dimensional  Model  of  the  Shell  Sub  Grid 
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4.  Cap  Sub  Grid 

The  cap  sub  grid  used  the  Lagrange  processor  and  was  filled  with  a  material 
model,  which  used  the  linear  equation  of  state,  the  Von  Mises  strength  model,  the  bulk 
strain  failure  model  and  the  erosion  model.  The  material  input  deck  can  be  found  in 
Appendix  D.  The  sub  grid  is  composed  of  %  of  a  solid  cylinder  with  an  outer  radius  of 
51.5mm  and  a  length  of  2mm  with  an  I  and  J  index  of  11  and  a  K  index  of  2  (75 
elements).  Figure  26  is  a  depiction  of  the  cap  sub  grid. 


Figure  26.  Three-Dimensional  Model  of  the  Shell  Sub  Grid 

The  cap  sub  grid  was  joined  to  the  shell  sub  grid  using  the  joined  Lagrange 
processor  coupling.  This  was  done  in  order  to  have  a  rigid  connection  between  the  steel 
sheet  casing  and  the  top  plate  of  the  mine. 
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V.  ANALYSIS  RESULTS 


All  of  the  analysis  results  for  this  thesis  were  obtained  using  a  Dell  Dimension 
XPS  T800r  computer  with  an  800  MHz  Intel  Pentium  Three  processor  and  640  mega 
bites  of  RAM. 

A.  FOUR-ROW  MINE  MODEL 

The  simulation  of  the  four-row  mine  model  required  a  run  time  of  approximately 
75  hours  and  36  minutes  and  a  total  of  10000  time  steps  were  computed,  with  a  final 
cycle  time  of  0. 1 59ms. 

After  0.074ms  (2900  cycles)  the  shell  sub  grid  and  the  wire  sub  grid  have 
expanded  in  the  radial  direction  with  the  shell  sub  grid  starting  to  break  apart  as  shown  in 
Figure  27.  With  Figure  28  showing  the  wire  and  shell  sub  grids  velocities  for  the  same 
cycle. 
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Figure  27.  Expansion  of  the  Four-Row  Mine  Model  at  0.074ms 
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Figure  28.  Four-Row  Mine  Model  Velocities  at  0.074ms 


By  0.159ms  (10000  cycles)  the  shell  sub  grid  has  significantly  deformed  and  is  in 
multiple  pieces  with  the  wire  sub  grid  radially  distributed  in  a  bunched  fashion.  The  final 
wire  and  shell  configuration  is  shown  in  Figure  29  with  the  final  velocities  shown  in 
Figure  30. 
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Figure  29.  Expansion  of  the  Four-Row  Mine  Model  at  0.159ms 


Velocity  vectors  (m/s) 
8.318E+005 


o  3 

o 

O  <3 


ozm-72 
Cycle  10000 
Time  1.593E-001  ms 
Units  mm,  mg,  ms 


o 

o  * 
3  O 


0.000E-H30C 


Figure  30.  Four-Row  Mine  Model  Velocities  at  0.159ms 
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The  target  point  locations  for  this  model  are  mapped  in  Figure  31.  At  these  target 
points  time  history  responses  were  computed.  The  pressure  applied  to  the  targets  and  the 
absolute  velocity  response  at  targets  1,  2,  3,  4,  5  and  at  targets  6,  7,  8,  9  are  plotted  and 
displayed  in  Figure  32.  The  peak  pressures  ranged  from  1  lxlCf  kPa  to  5x10s  kPa  with 
the  lower  peak  pressures  occurring  at  the  outer  most  targets.  The  pressure  increased  to 
the  peak  pressure  almost  instantaneously  and  then  decayed  exponentially  as  time 
increased.  This  is  what  should  be  expected  since  an  explosive  shock  wave  was  used  in 
modeling  the  mine  explosion.  The  peak  velocities  ranged  from  200  m/s  to  450  m/s.  The 
lower  velocities  tend  to  be  the  inner  most  fragments  while  the  outer  fragments  had  the 
higher  velocities.  The  oscillation  in  the  pressure  and  the  velocity  plots  appears  to  have 
been  caused  by  the  fragments  coming  in  contact  with  the  symmetry  boundary  along  the  z- 
axis. 


Figure  31.  Four-Row  Mine  Model  Target  Point  Locations 
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Figure  32.  Four-Row  Mine  Model  Pressure  and  Velocity  Time  History  Plots 


B.  THREE-ROW  MINE  MODEL 

The  simulation  of  the  three-row  mine  model  required  a  run  time  of  approximately 
8  hours  and  41  minutes  and  a  total  of  6000  time  steps  were  computed,  with  a  final  cycle 
time  of  0.201ms. 

After  0.103ms  (3000  cycles)  the  shell  sub  grid  and  the  wire  sub  grid  have 
expanded  in  the  radial  direction  as  shown  in  Figure  33.  With  Figure  34  showing  the  wire 
and  shell  sub  grids  velocities  for  the  same  cycle. 
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Figure  33.  Expansion  of  the  Three-Row  Mine  Model  at  0.103ms 


Cycle  3000 

Time  1.033E-001  ms 

Units  mm,  mg,  ms 


Velocity  vectors  (m/s) 
6. 082  E +002 


5.474E+00: 

4. 866  E +00: 

4.257E+00: 

3.649E+00: 

3.041  e+oo: 

2.433E+00: 

1 .825E+00: 

1.216E+00: 

6  082E+00 
O.OOOE+Ooj 


l 


Figure  34.  Three-Row  Mine  Model  Velocities  at  0.103ms 

36 


By  0.201ms  (6000  cycles)  the  shell  sub  grid  has  started  to  break  apart  and  is  in 
several  pieces.  The  wire  sub  grid  is  radially  distributed  in  a  fairly  symmetric  manner. 
The  final  wire  and  shell  configuration  is  shown  in  Figure  35  with  the  final  velocities 
shown  in  Figure  36. 


Figure  35.  Expansion  of  the  Three-Row  Mine  Model  at  0.2ms 
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Figure  36.  Three-Row  Mine  Model  Velocities  at  0.2ms 


The  target  point  locations  for  this  model  are  mapped  in  Figure  37.  At  these  target 
points  time  history  responses  are  computed.  The  pressure  applied  to  the  targets  and  the 
absolute  velocity  response  at  targets  6,  11,  16  and  at  targets  8,  13,  22,  23  are  plotted  and 
displayed  in  Figure  38.  The  peak  pressures  ranged  from  12xl05  kPa  to  7xl05  kPa  with 
the  lower  pressures  occurring  at  the  outer  fragments.  In  general  the  pressure  increases  to 
the  peak  pressure  almost  instantaneously  and  then  decays  exponentially  as  time  increase. 
This  is  what  was  expected  since  an  explosive  shock  wave  was  used  in  modeling  the  mine 
explosion.  The  peak  velocities  ranged  from  225  m/s  to  300  m/s.  The  lower  velocities 
tend  to  be  the  inner  most  fragments  while  the  outer  fragments  have  higher  velocities.  The 
sudden  drops  and  increases  in  the  fragments  velocities  seams  to  correspond  to  the 
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Pressure  (kPa) 


fragments  hitting  each  other.  These  changes  in  velocity  also  correspond  to  the  spikes  in 
the  pressure  plots. 


Figure  37.  Three-Row  Mine  Model  Target  Point  Locations 


Figure  38.  Three-Row  Mine  Model  Pressure  and  Velocity  Time  History  Plots 
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C.  TWO-ROW  MINE  MODEL 

The  simulation  of  the  two-row  mine  model  required  a  run  time  of  approximately  2 
hours  and  32  minutes  and  a  total  of  7000  time  steps  were  computed,  with  a  final  cycle 
time  of  0.216ms. 

After  0.12ms  (3500  cycles)  the  shell  sub  grid  and  the  wire  sub  grid  have  expanded 
radially  with  the  shell  sub  grid  starting  to  break  apart  as  shown  in  Figure  39  with  Figure 
40  showing  the  wire  and  shell  sub  grid  velocities  for  the  same  cycle. 
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Figure  39.  Expansion  of  the  Two-Row  Mine  Model  at  0.12ms 
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Figure  40.  Two-Row  Mine  Model  Velocities  at  0.12ms 


By  0.216ms  (7000  cycles)  the  shell  is  completely  fragmented  and  is  intermingled 
with  the  wire  fragments.  The  fragments  are  fairly  evenly  distributed  with  several  small 
groups  of  fragments.  The  final  wire  and  shell  configuration  is  shown  in  Figure  41  with 
the  final  velocities  shown  in  Figure  42. 
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Figure  41.  Expansion  of  the  Two-Row  Mine  Model  at  0.216ms 
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Figure  42.  Two-Row  Mine  Model  Velocities  at  0.216ms 
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The  target  point  locations  for  this  model  are  mapped  in  Figure  43.  At  these  target 
points  time  history  responses  are  computed.  The  pressure  applied  to  the  targets  and  the 
absolute  velocity  response  at  targets  1,  2,  3  and  at  targets  8,  9,  10  are  plotted  and 
displayed  in  Figure  44.  The  peak  pressure  ranged  from  8xl05  kPa  to  5xl05  kPa.  The 
pressure  increases  to  the  peak  pressure  almost  instantaneously  and  then  decays 
exponentially  as  time  increase.  This  is  what  should  be  expected  since  an  explosive  shock 
wave  was  used  in  modeling  the  mine  explosion.  The  peak  velocities  ranged  from  225 
m/s  to  250  m/s.  The  velocities  of  the  targets  for  this  model  were  for  the  most  unifonn 
with  them  within  25  m/s  of  each  other.  In  addition,  discontinuities  in  the  velocity  and 
pressure  plots  occurred.  These  appear  to  have  been  caused  by  the  fragments  colliding 
with  each  other  or  the  shell  sub  grid. 


Figure  43.  Two-Row  Mine  Model  Target  Point  Locations 
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Figure  44.  Two-Row  Mine  Model  Pressure  and  Velocity  Time  History  Plots 

D.  THREE-DIMENTIONAL  MINE  MODEL 

The  simulation  of  the  three-dimensional  mine  model  required  a  run  time  of 
approximately  89  hours  and  32  minutes  and  a  total  of  4500  time  steps  were  computed, 
with  a  final  cycle  time  of  0. 149ms. 

After  0.7ms  (2200  cycles)  the  shell  sub  grid  and  the  wire  sub  grid  have  expanded 
radially  with  the  cap  sub  grid  separated  from  the  shell  sub  gird  as  shown  in  Figure  45. 
The  velocities  of  the  wire,  shell  and  cap  sub  grids  for  the  same  cycle  are  shown  in  Figure 
46. 
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Figure  45.  Expansion  of  the  Three-Dimensional  Mine  Model  at  0.07ms 
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Figure  46.  Three-Dimensional  Mine  Model  Velocities  at  0.07ms 
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By  0.149ms  (45000  cycles)  the  shell  sub  grid  has  broken  and  the  cap  sub  gird  is 
significantly  deformed.  The  final  wire,  shell  and  cap  configurations  are  shown  in  Figure 
47  with  the  final  velocities  shown  in  Figure  48. 


Units  mm,  mg,  ms 

Figure  47.  Expansion  of  the  Three-Dimensional  Mine  Model  at  0.149ms 
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Figure  48.  Three-Dimensional  Mine  Model  Velocities  at  0.149ms 


The  target  point  locations  for  this  model  are  mapped  in  Figure  49.  At  these  target 
points  time  history  responses  are  computed.  The  pressure  applied  to  the  targets  and  the 
absolute  velocity  response  at  targets  17,  18,  19  and  at  targets  24,  25,  26  are  plotted  and 
displayed  in  Figure  50.  These  targets  are  the  fragments  closest  to  the  center  of  the  mine. 
The  peak  pressure  ranged  from  7.5xl05  kPa  to  5x10s  kPa.  In  general  these  plots  compare 
well  with  that  of  the  two-row  mine  model.  This  is  expected  since  this  model  was  an 
extension  of  that  model.  The  peak  velocities  ranged  from  190  m/s  to  240  m/s.  The  lower 
velocities  tend  to  be  the  inner  most  fragments  while  the  outer  fragments  have  higher 
velocities.  In  general  these  plots  compare  well  with  that  of  the  two-row  mine  model. 
This  is  expected  since  this  model  was  an  extension  of  that  model. 
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Figure  49.  Three-Dimensional  Mine  Model  Target  Points 
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Figure  50.  Three-Dimensional  Mine  Model  Pressure  and  Velocity  Time  History  Plots 
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The  next  group  of  target  points  plotted  are  for  the  targets  34,  35,  36  and  41,  42  , 
43.  These  targets  are  the  fragments,  which  are  Vi  of  the  way  between  the  top  and  the 
middle  of  the  mine.  The  peak  pressure  ranged  from  7.5xl05  kPa  to  4.9xl05  kPa  with  the 
peak  velocities  ranging  from  175  m/s  to  225  m/s.  These  time  history  plots  can  be  found 
in  Figure  51. 


Figure  51.  Three-Dimensional  Mine  Model  Pressure  and  Velocity  Time  History  Plots 


The  final  group  of  target  points  plotted  are  for  the  targets  51,  52,  53  and  58,  59, 
60.  These  targets  are  the  fragments,  which  are  at  the  top  of  the  mine  near  the  cap  sub 
grid.  The  peak  pressure  ranged  from  6X103  kPa  to  4x10^  kPa  with  the  peak  velocities 
ranging  from  150  m/s  to  175  m/s.  These  time  history  plots  can  be  found  in  Figure  52. 
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The  reason  for  this  decrease  in  pressure  and  velocity  is  due  to  the  fact  that  mine  model  is 
not  constrained  in  the  vertical  direction.  A  portion  of  the  pressure  and  energy  is  used  to 
propel  the  cap  sub  grid  upward,  reducing  the  effect  on  the  fragments  near  the  cap  sub 
grid. 


Figure  53.  Three-Dimensional  Mine  Model  Pressure  and  Velocity  Time  History  Plots 


An  estimation  of  the  probability  of  kill  was  computed  using  the  three-dimensional 
models  data  to  compute  the  results.  The  probability  of  kill  was  computed  using  the  Mat 
Lab  program  provided  in  Appendix  E.  This  program  was  developed  based  on  the 
equations  and  assumptions  outlined  in  Chapter  III  of  this  report.  In  Figure  54  the 
probability  of  kill  is  plotted  with  respect  to  distance  for  a  soft  target.  Each  line  in  the 
figure  represents  the  targets  width,  which  range  from  0.05  m  to  1  m.  In  Figure  55  the 
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same  type  of  plot  is  produced  for  the  hard  target.  When  comparing  the  two  plots  it  is 
apparent  that  for  the  harder  target  the  probability  of  kill  decays  at  a  considerably 
increased  rate  than  that  of  the  soft  target.  This  is  due  to  the  fact  that  the  harder  target  was 
able  to  absorb  more  of  the  fragments  kinetic  energy. 


Figure  54.  Estimated  Probability  of  Kill  for  a  Soft  Target 


Figure  55.  Estimated  Probability  of  Kill  for  a  Hard  Target 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  coupled  Lagrangian  and  Eulerian  approach  was  used  to  model  the  detonation 
and  fragmentation  phenomena  in  both  two  and  three-dimensions.  The  analysis  results 
from  these  models  demonstrated  that  the  maximum  fragment  velocity  is  extremely  high, 
releasing  high  kinetic  energy.  In  addition  the  probability  of  kill  for  soft  and  hard  targets 
was  estimated  as  a  function  of  the  stand  off  distance.  This  revealed  that  the  rate  of  decay 
of  the  probability  of  kill  for  the  hard  target  was  considerably  faster  than  that  of  the  soft 
target.  The  coupled  Lagrangian  and  Eulerian  approach  may  be  adequate  for  the  modeling 
to  the  fragmentation  phenomena. 

It  is  recommended  that  additional  studies  be  conducted  into  the  detonation  and 
fragmentation  phenomena.  These  studies  include  but  are  not  limited  to  an  investigation 
into  what  effect  the  resin  that  holds  the  fragments  in  place  has  on  the  fragmentation 
effect.  In  addition  the  development  of  equations  of  state  and  material  models  that  take 
thennal  effects  into  account  could  be  beneficial  in  modeling  this  phenomena. 
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APPENDIX  A.  MATIERIAL  INPUT  DECK  FOR  TNT  SUB  GRID 


MATERIAL  NAME:  TNT 
EQUATION  OF  STATE:  JWL  (Explosive) 

Reference  density  (g/cm3):  1.63000E+00 
Parameter  Ci  (kPa):  3.73770E+08 
Parameter  C7  (kPa):  3.74710E+06 
Parameter  Rl:  4.15000E+00 
Parameter  R2:  9.00000E-01 
Parameter  W:  3 . 5  OOOOE-0 1 
C-J  Detonation  velocity  (m/s):  6.93000E+03 
C-J  Energy  /  unit  volume  (kJ/m3):  6.00000E+06 
C-J  Pressure  (kPa):  2.10000E+07 
Bum  on  compression  fraction:  0.00000E+00 
Pre-bum  bulk  modulus  (kPa):  0.00000E+00 
Adiabatic  constant  (kPa):  0.00000E+00 
STRENGTH  MODEL:  None  (Hydro) 

FAILURE  MODEL:  None 
EROSION  MODEL:  None 
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APPENDIX  B.  MATERIAL  INPUT  DECK  FOR  AIR  SUB  GRID 


MATERIAL  NAME:  AIR 
EQUATION  OF  STATE:  Ideal  Gas 
Gamma:  1 .40000E+00 
Reference  density  (g/cm3):  1.22500E-03 
Adiabatic  Constant  (g/cm3):  0.00000E+00 
Pressure  shift  (kPa):  0.00000E+00 
Reference  Temperature  (K):  2.88200E+02 
Specific  Heat  (C.V.)  (J/kgK):  7.17300E+02 
STRENGTH  MODEL:  None  (Hydro) 

FAILURE  MODEL:  None 
EROSION  MODEL:  None 
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APPENDIX  C.  MATERIAL  INPUT  DECK  FOR  WIRE  SUB  GRID 


MATERIAL  NAME:  WIRE 
EQUATION  OF  STATE:  Linear 

Reference  density  (g/cm3):  7.90000E+00 
Bulk  Modulus  (kPa):  2.00000E+08 
Reference  Temperature  (K):  0.00000E+00 
Specific  Heat  (C.V.)  (J/kgK):  0.00000E+00 
STRENGTH  MODEL:  Elastic 

Shear  Modulus  (kPa):  9.00000E+07 
FAILURE  MODEL:  None 
EROSION  MODEL:  None 
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APPENDIX  D.  MATERIAL  INPUT  DECK  FOR  SHELL  AND  CAP 

SUB  GRIDS 


MATERIAL  NAME:  STEEL 
EQUATION  OF  STATE:  Linear 

Reference  density  (g/cm3):  7.90000E+00 
Bulk  Modulus  (kPa):  2.00000E+08 
Reference  Temperature  (K):  0.00000E+00 
Specific  Heat  (C.V.)  (J/kgK):  0.00000E+00 
STRENGTH  MODEL:  Von  Mises 

Shear  Modulus  (kPa):  9.00000E+07 
Yield  Stress  (kPa):  2.00000E+05 
FAILURE  MODEL:  Bulk  Strain 

Ultimate  Strain:  1.50000E-01 
Crack  Softening,  Gf  (J/m2):  0.00000E+00 
or,  Kc  (mN/m3/2):  0.00000E+00 
EROSION  MODEL:  Inst.  Geo.  Strain 
Erosion  Strain:  2.00000E+00 
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APPENDIX  E.  PROBABILITY  OF  KILL  PROGRAM 


The  following  program  was  developed  to  estimate  the  probability  of  kill  for  the 
mine  models  developed  in  this  study.  The  program  uses  the  equations  and  assumptions 
outlined  in  Chapter  III  of  this  report. 

load  fragdata 


%  Constants 
rho  =  1.2; 

Cd  =  2.1; 

A  =  0.000025; 
m  =  0.0031; 

S  =  [1,  5,  10,  15,25]; 

N  =  [1224]; 

TA  =  [1,  .75,  .5,  .25,  .1,  .05]; 


%  Air  Density  at  Standard  Conditions 
%  Drag  Coefficient  of  a  Square  Cylinder 
%  Cross-Sectional  Area  of  a  Fragment 
%  Mass  of  a  Fragment 
%  Distance  from  Mine 
%  Number  of  Fragments  in  3D  Mine  Model 
%  Width  of  Target 


%  Determine  the  average  initial  velocity  of  the  fragments. 
Vor3d  =  mean(R3D2VEL(91,:)); 


%  Determine  the  velocity  of  fragments  at  specified  distances  considering  drag, 
for  n  =  1:5 

Vsr3d(n,:)  =  Vor3d  .*  exp(-(rho  *  A  *  S(n))/(2  *  m));  %  Equation  3.2 

FragD(n)  =N/(2  *  pi  *  S(n)); 
end 


%  Determine  kinetic  energy  of  projectiles  based  on  velocity  at  specified  distances. 
KEr3d  =  (Vsr3d.A2)  .*  (1/2)  *  m;  %  Equation  3.1 


%  Determine  Probability  of  Kill  for  Targets 
Pkhl  =  (KEr3d  ./  100)  .*  0.1; 

Pkh2  =  (KEr3d  ./  4000)  .*  0.1; 
for  i  =  1:5 
for  j  =  1:6 

Nd(j,i)  =  FragD(i)  *  TA(j);  %  Equation  3.5 

if  Nd(j,i)  >  1 

PKl(j,i)  =  1  -(l-Pkhl(i))ANd(j,i);  %  Equation  3.3 
PK2(j,i)  =  1  -(l-Pkh2(i))ANd(j,i); 
else 

PKl(j,i)  =Nd(j,i)  *  Pkhl(i);  %  Equation  3.4 

PK2(j,i)  =  Nd(j,i)  *  Pkh2(i); 
end 
end 
end 

figure  (1) 
plot(S',PKl) 

legend  (Tm1, '0.75m', '0.5m', '0.25m', '0.  lm', '0.05m') 

figure  (2) 

plot(S',PK2) 

legend  ('lm', '0.75m', '0.5m', '0.25m', '0.  lm', '0.05m') 


%  Soft  Target 
%  Flard  Target 
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